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for a global bifurcation 
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We report an experimental evidence ol a global bilurcation on a highly turbulent von Karman flow. 
The mean flow presents multiple solutions: the canonical symmetric solution becomes marginally 
unstable towards a flow which breaks the basic symmetry of the driving apparatus even at very 
large Reynolds number. The global bifurcation between these states is highly subcritical and the 
system thus keeps a memory of its history. The transition recalls low-dimension dynamical systems 
transitions and exhibits a very peculiar statistics. We discuss the role of turbulence in two ways: the 
multiplicity of hydrodynamical solutions and the effect of fluctuations on the nature of transitions. 

PACS numbers: 05.45.-a, 47.20.-k, 47.27.Sd 



Non- linear systems generally present multiple solutions 
and various transitions between them. Moreover, stabil- 
ity and transitions are influenced by the presence of noise 
and/or fluctuations. In the field of turbulence, the ques- 
tion of multistability of turbulent flows, for example in 
tornadoes [10] , delta wing flow Q , wakes Q , and vortex 
breakdown |5|], remains open and unsolved. While mul- 
tiple analytical or numerical solutions are often encoun- 
tered for the Navier-Stokes equation at even moderate 
Reynolds number (e.g., for swirling flows 0, HI IE 01 E3] ) j 
these solutions are generally neither experimentally rel- 
evant, nor stable at very high Reynolds number. Fur- 
thermore, turbulent flows at very high Reynolds number 
are generally expected to statistically respect the basic 
symmetries of their driving apparatus. Indeed, even if 
bifurcations and symmetry breaks occur on the way to 
turbulence, the fully developed turbulent state is known 
to restore the broken symmetries, in the limit of infi- 
nite Reynolds number and far from boundaries [9|- In 
this Letter, we experimentally study the multistability 
of a turbulent von Karman flow between two counter- 
rotating disks in a finite vessel at very high Reynolds 
number. This system undergoes a subcritical global bi- 
furcation between turbulent states characterized by mean 
flows of different topology and symmetry. These turbu- 
lent states coexist at high Reynolds number and can be 
"prepared" specifically, i.e., they keep a memory of the 
system history. Since these states are highly fluctuat- 
ing turbulent states, we also address the question of the 
role of the fluctuations for such a transition. Actually, 
the effect of an external noise on an existing transition is 
well documented ^3 , but the global bifurcation reported 
here does only take place over an already fluctuating tur- 
bulent regime. Do fluctuations trigger the bifurcation as 
multiplicative noise do for nonlinear oscillators ^lj an d 
turbulent a-effect do for dynamo action ^3] ? 

Experimental setup. We call von Karman type flow 
the flow generated between two coaxial counter-rotating 
impellers in a cylindrical vessel. The cylinder radius and 
height are respectively R = 100 mm and H c = 500 mm. 



We use bladed disks to ensure inertial stirring. Most of 
the inertially driven von Karman setups studied in the 
past dealt with straight blades 0, 0j. In this Letter, 
the impellers consist of 185 mm diameter disks each fitted 
with 16 curved blades — curvature radius 50 mm, height 
20 mm (Fig. ^ . The distance between the inner faces of 
the disks is H = 180 mm which defines a working space 
for the flow of aspect ratio H/R = 1.8. With curved 
blades, the directions of rotation are no longer equivalent. 
We rotate the impellers clockwise (with the concave face 
of the blades). Four baffles (10 x 10 x 125 mm) can be 
added along the cylinder wall. 

The impellers are driven by two independent brush- 
less 1.8 kW motors, with a speed servo loop control. The 
motor rotation frequencies /i, /a can be varied indepen- 
dently in the range — 15 Hz. An experiment is thus char- 
acterized by two numbers: / = y/ (/j + /f )/2 measuring 
the intensity of the forcing and 9 = (fa — /i)/(/i + /a) 
measuring the speed dissymmetry (—1 < 9 < 1). For ex- 
act counter rotation /i = /2 = / and 9 = 0. The speed 
servo loop control ensures a precision of 0.5% on /, and 
an absolute precision of ±0.002 on 9 for small values. 

The working fluid is water. Copper cooling coils be- 
hind the impellers and a thermoregulated bath ensure 
a thermal regulation with a precision of 1°C. Velocity 
fields are measured by Laser Doppler Velocimetry (LDV) . 
Torques are measured as an image of the current con- 
sumption in the motors given by the servo drives and 
have been calibrated by calorimetry. The analog signal 
is low-pass filtered at 10 Hz. For a typical frequency 
/ = 4 Hz at 35°C, the integral Reynolds number is 
Re = 2nfR 2 v~ 1 ~ 3. 10 5 and the velocity fluctuation 




FIG. 1: Sketch of the experimental setup and of the impellers 
blades profile. The arrow indicates the positive rotation sense. 
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level is of order 30%: the flow is highly turbulent. 

The von Karman flow phenomenology is the follow- 
ing. Each impeller acts as a centrifugal pump: the fluid 
rotates with the impeller and is expelled radially, ft 
is pumped in the center of the impeller. In the exact 
counter-rotating regime, the flow is divided into two toric 
cells separated by an azimuthal shear layer. The problem 
(equation and boundary conditions) is invariant under ro- 
tations of it (72-Tr) around any radial axis passing through 
the center of the cylinder. The velocity field is expected 
T^-invariant. 

A "statistical" symmetry breaking. In our high 
Reynolds number regime, the flow is highly turbulent. 
For instance the rms value of the velocity is of the same 
order of magnitude as the mean value. In Fig. (left), 
we present a map of the mean part of the exact counter- 
rotation flow measured by LDV. Two cells are observed, 
the flow is TZtt -invariant: the symmetries are statistically 
restored [jj . The mean angular momentum of the fluid is 
equal to zero: the two impellers produce the same mean 
torque to maintain the flow. This situation is well-known 
and documented. We label this symmetric state (s). 

However, with our curved blades, we observe for small 
9 a global bifurcation of the flow after a certain time tu / : 
both mean velocity field and torques display dramatic 
changes (Fig. OJ. The two torques are suddenly 4 times 
larger, and are no longer equal. The mean flow exhibits 
only one cell (Fig. right). In the bulk, the fluid is 
pumped toward impeller 1 without rotation. Then the 
fluid is expelled radially and starts spiraling along the 
cylinder until it meets impeller 2 which rotates in the 
opposite direction. It is abruptly stopped and reinjected 
near the axis. We label this state (bi). A third state (62) 
is deduced from (61 ) by exchanging the roles of impellers 
1 and 2. In bifurcated states (61) or (62), the fluid is 
globally in rotation: the mean angular momentum is not 
zero. 

Finally, three states are observed: the canonical TZ^- 
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FIG. 2: Dimensionless mean velocity field measured at 9 — 
by LDV over 120 integral turn-over time by grid point to 
ensure good convergence; f — 2 Hz (Re = 1.5 10 ). Left: 
symmetric state (s). Right: bifurcated state (61). Space co- 
ordinates in units of R. Gray code stands for azimuthal ve- 
locity. Isolines are distant of 0.2 and the gray code saturates 
in the right map. Bold lines indicate level zero. 



invariant — in a statistical sense — state (s) and two bifur- 
cated states which break the TZ n symmetry at 9 = 0, but 
are images one of the other by TZ^ . We detail in the next 
section the transitions between these different states. 

Hysteresis loops. The difference between the two 
torques characterizes the different states. We have 
checked that, as expected for so high a Reynolds num- 
ber 0, the torque T given by one impeller for a given 
(/, 9) does not depend on Re and scales as: T(f, 9) = 
K, p {9) pR 5 (2tt/) 2 14], with p the fluid density and K p 
a dimensionless power coefficient. 

In Fig. 01 we plot the dimensionless difference AK p 
between the two torques vs. 9 for several configurations. 
For straight blades, we observe a continuous curve from 
9 = — 1 to 6 = 1 (Fig. HJi) with two transitions between 
one- and two-cells flows at 9 = ±0.13. For impellers with 
curved blades and no baffles on the cylinder wall, we ob- 
serve the three states in Fig. 0Jx For 9 = 0, we recognize 
state (s) (AK p = 0), and both bifurcated states (61) and 
(62). State (s) branch is almost reduced to one point and 
can only be reached by starting the two motors simul- 
taneously. Its stability is discussed in the next section. 
The bifurcated state (61) lies on a branch coming contin- 
uously from 9 = — 1 (/2 = 0). Starting from 9 = — 1 and 
increasing 9, we stay on the (61) branch even for 9 > 0: 
impeller 1 keeps rotating and pumping the fluid although 
its rotation rate is weaker than impeller 2 rotation rate. 
For 9 ~ 0.16 there is a transition from (bi) to (62): the 
fluid abruptly changes its sense of rotation. There is a 
large hysteretic cycle. Note that it is impossible to reach 
the symmetric state (s) by this way. The global quanti- 
ties of this highly turbulent flow keep memory of the way 
the system has been started from rest. An intermedi- 
ate situation is reached with the same curved blades and 
baffles on the wall (Fig. |3J;). Baffles break the spiraling 
flow along the wall of the cylinder, which is a major fea- 
ture of the bifurcated state velocity field. The hysteretic 
cycle splits into two classical first order cycles: the cen- 
tral symmetric state becomes stable and can be obtained 
from any initial condition. 

Stability of the central branch (s). We focus now on 
the transition from symmetric state to bifurcated state 
for curved blades without baffles. As mentioned before, 
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FIG. 3: Time serie of dimensionless torque showing the bifur- 
cation (s) -> (62), for 6 = 0.0204, / = 4.08 Hz. Left: torque 
on impeller 1. Right: torque on impeller 2. The bifurcation 
time is the time when the torque on impeller 1 reaches 140% 
of the mean value for the symmetric state (s). 
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the central branch is very small and, for a given (/, 9), 
the transition occurs after a certain time tbif which ex- 
hibits a complex statistics. 

So we performed the following experiments: starting 
from rest, we simultaneously start both motors to a given 
(/, 9) with a short ramp (typically 1 s) and record the 
torques. Few seconds after the instant tut when bifur- 
cation occurs, we stop the motors, wait a minute and 
run again. We perform typically 500 runs to get the 
distribution of bifurcation times. The cumulative dis- 
tribution function (CDF) for tuf (Fig. EJ) shows ex- 
ponential behavior for the probability of staying in the 
symmetric state a time greater than t: P{tbif > t) = 
A exp[—(t — to)/r], to is characteristic of the transition 
duration (to-/ ~ 5). Thus, we obtain a characteristic bi- 
furcation time t(/, 9) by non-linear fitting of the CDF. 
We performed the experiment for three values of /. The 
results are shown in Fig. H3 in log-log scale. There is no 
noticeable dependence on / and t behaves as | 9 |~ 6 . 
So, as 9 tends to zero, r diverges very fast to infinity: 
the central point is marginally stable. The physical phe- 
nomenon at the origin of such an exponent remains to be 
understood. 
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FIG. 4: Dimensionless torque difference AK P vs. 8 for Re 
in the range 2 — 8 x 10 . Straight blades (a) exhibit con- 
tinuous transition from 1-cell flow to 2-cells flow for 8 = 
±0.13 (vertical lines). Curved blades without (b) or with 
(c) baffles along cylinder wall show subcritical transitions be- 
tween symmetric/2-cells (s)-(D) and bifurcated/ 1-cell (6i)- 
(6 2 ) states (O). 
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FIG. 5: Cumulative density function of bifurcation times for 
three different 8 at / = 4.16 Hz. Dotted line: non-linear 
exponential fit. 



Discussion. The experiment presented here opens 
mainly two problems: (i) the existence and the nature 
of multiple regimes for this turbulent flow, and (ii) the 
role of the noise or the fluctuations in some transitions 
between these flow regimes, i.e., the stability problem for 
the two-cell (s) branch. 

We first try to explain the existence of multiple stable 
regimes by hydrodynamical basic arguments. The von 
Karman (VK) class of Navier-Stokes solutions in semi- 
infinite space with one or two infinite rotating disks for 
end-conditions has been extensively studied since 1921 
Q> ^3 • Experiments are necessarily limited in diame- 
ter and do not strictly belong to the same class. However 
the approximation is very commonly made at least for 
small H/R. In practice, in our system, and in the spirit 
of Batchelor and Stewartson [tJ , we construct finite- 
aspect-ratio solutions of our experimental VK problem at 
high Re (Fig.[5J with (i) any typical truncated Batchelor 
6] solution for < r < R/2 together with (ii) some recir- 
culation flow in rotation in R/2 < r < R and (iii) a thin 
boundary layer near the outer cylinder which matches 
this rotation [l^. The two-cell mean flow (s) is sim- 
ply described in the laboratory frame by two rotating 
regions inertially driven by the blades and separated by 
a shear layer near mid- height. Both disks centrifugally 
expel the fluid. Let's now consider one-cell flows (b\) and 
(62). Since one disk expels the fluid and the other rein- 
jects it to the center, these flows resemble the corotating 
(fi-fz < 0) regime solutions 0,0 characterized by: uni- 
form rotation of the bulk; a boundary layer on each disk; 
pumping from one disk to the other and recirculation 
at infinity. This solution has no shear layer. Let's note 
that mean bulk rotation — for r < R/2 — is close to zero 
(Fig. 13 right). In conclusion, we can make the assump- 
tion that flows (61) and (62) are equivalent to corotating 
flows observed in two oppositely rotating frames of fre- 
quencies +f r and —f r with |/ r | > max(|/i |, |/ 2 |). This 
is well consistent with the fact that the single cell flows 
(61) and (62) exhibit global rotation in the outer shell 
;$ t < R. The stability of such solution is clearly 
enhanced by the concave curved blades that enforce ro- 
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FIG. 6: T.f vs. 8 for / = 4.16 Rz/Re = 3.3 10 5 (O), / = 6 
Hz/ Re = 4.7 10 5 (□) and / = 10 Hz/ Re = 7.9 10 5 (0), fitted 
by a —6 slope power law. 
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tation of the fluid near the outer cylinder. 

Let's now consider how these three solution branches 
exchange their stability. First note that the bifurcation 
diagrams respect the IZ^ symmetry: 9 — ► —9; AK p — ► 
—AK p . The straight blades diagram (Fig. is con- 
tinuous: from left to right two second order transitions 
(bi)<->(s) and (s)<->(&2) are observed as in small H/R 
systems 0. On the contrary, the curved blades dia- 
gram (Fig. Et>) is strongly hysteretic. Addition of baffles 
(Fig. Efc) allows to remove a degeneracy: baffles drag dis- 
turb the outer cylinder boundary layer flow, thus lower- 
ing the relative stability of one-cell flows with respect to 
the two-cell flow. The large hysteresis cycle is split into 
two classical first-order bifurcations. This singular cycle 
can thus be viewed as the result of the collapse or collision 
of two first-order cycles. Similar cycles are encountered 
in conical 0] and delta- wing flows Q . The memory ef- 
fect — if the system is currently on (s), both driving fre- 
quencies must have been increased in parallel — is thus 
essentially a consequence of the cycle structure. 

In order to test the effect of turbulence on the stabil- 
ity of the observed flows, we lowered Reynolds number 
down to laminar using water/glycerol mixtures. While 
Re < 1000, no multiplicity is observed: the bifurca- 
tion diagram is similar to the straight-blade diagram of 
Fig-UK- The cycle appears for Re between 1000 and 3000. 
The study is in progress and will be reported elsewhere. 
The high Reynolds behavior reported in this Letter is well 
established once Re > 5000. Thus, multiplicity appears 
with turbulence and does not with laminar (Re < 110) 
nor chaotic (Re < 1000) flows. A possible explanation 
for the multiplicity could thus be the evolution of the 
outer cylinder boundary layer with Re. 

Besides, the statistical nature of the transitions them- 
selves is probably related to turbulent fluctuations. Let's 
first notice that the bifurcation studied here corresponds 
to exchange of stability between mean flows, these mean 
states being never realized at any given time. Is the bi- 
furcation formalism exactly valid for our mean flows ? 
On some aspects, our system behaves as a low-dimension 
dynamical system, as in the turbulent spiral transition 
observed in wide-aspect-ratio Taylor-Couette flow 0] or 
in the noise-induced Hopf bifurcation for a Duffing os- 
cillator with multiplicative white noise |ll| . However, 
suppose a non-linear amplitude equation could correctly 
describe the shape bifurcation diagram, it would proba- 
bly not be able to catch the statistics of the transition 
from the two-cell state (s) to a one-cell state (b\) or (62)- 
This transition shows a very peculiar statistics, with a 
very high critical exponent 6 (Fig. |BJ). It also strictly 
respects a forbidden-transition rule: the horizontal axis 
of the bifurcation diagram is never crossed, i.e., (s)— +(62) 
[resp. (s)— ►(&].)] is forbidden for 9 < [resp. 6 > 0]. This 
observational fact could by itself justify the stability of 
the central point 6 = 0, which has to respect both rules. 
Furthermore, the non-crossing of the axis could be the 



signature of multiplicative noise as suggested to account 
for small-scale turbulence [20| . 

The global bifurcation reported in this Letter presents 
a very unusual bifurcation diagram. Some features about 
the multistability have been searched among the mechan- 
ics of high- Reynolds- number flows, while some other sim- 
ply involve the theory of non-linear bifurcations, possi- 
bly in the presence of noise. Among the transitions, the 
two-cell — > one-cell stability exchange plays a remarkable 
role, in presenting an original statistics of transition and 
putting the flow definitively in a state which breaks the 
TZ„ symmetry of the system and does not allow the flow 
to restore statistically this symmetry when Re — > 00. 
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